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ABSTRACT 

A number of relaxed, cool-core galaxy clusters exhibit diffuse, steep-spectrum radio sources in 
their central regions, known as radio mini-halos. It has been proposed that the relativistic electrons 
responsible for the emission have been reaccelerated by turbulence generated by the sloshing of the 
cool core gas. We present a high-resolution MHD simulation of gas sloshing in a galaxy cluster coupled 
with subgrid simulations of relativistic electron acceleration to test this hypothesis. Our simulation 
shows that the sloshing motions generate turbulence on the order of Sv ~ 50-200 km s~^ on spatial 
scales of ~50-100 kpc and below in the cool core region within the envelope of the sloshing cold fronts, 
whereas outside the cold fronts, there is negligible turbulence. This turbulence is potentially strong 
enough to reaccelerate relativistic electron seeds (with initial 7 ~ 100 — 500) to 7 ~ 10^ via damping 
of magnetosonic waves and non-resonant compression. The seed electrons could remain in the cluster 
from, e.g., past AGN activity. In combination with the magnetic field amplification in the core, these 
electrons then produce diffuse radio synchrotron emission that is coincident with the region bounded 
by the sloshing cold fronts, as indeed observed in X-rays and the radio. The result holds for different 
initial spatial distributions of preexisting relativistic electrons. The power and the steep spectral index 
(a ~ 1 — 2) of the resulting radio emission are consistent with observations of minihalos, though the 
theoretical uncertainties of the acceleration mechanisms are high. We also produce simulated maps 
of inverse- Compton hard X-ray emission from the same population of relativistic electrons. 
Subject headings: galaxies: clusters: general. X-rays: galaxies: clusters, turbulence, MHD, radio 
continuum: galaxies 



1. INTRODUCTION 



A number of relaxed, cool-core clusters are hosts 
to faint, diffuse radio emission with a radius com- 
parable to the size of the cooling region (r ^ 100- 
300 kpc) and a steep spectrum {a > 1; S^ oc v~^). 
These sources, called mini-halos, are relatively rare, 
with currently only around 10 clusters wit h confirmed 
detec tions. Examples incl ude Perseus (iBurns et al.l 
[T99I [sKbring 1993|) A2029 (iGovoni et aT]|2QQ9|) , OphT 
uchus^Govoni etldl '20 091: iMur^iaet al.l [20To[ ). RXC 
J1504T^^248_(Giacintucci^elalJ|20nl), and RXJ 1347- 
1145 (|Gitti et al.ii2007 ). to name a few. Questions still 
remain about the physical properties and the origin of 
these sources. 

Though clusters hosting mini-halos have central active 
galaxies, they are not sufficient by themselves to power 
the diffuse radio emission. The radiative timescale of 
the electrons at the required energies for the observed 
emission (~ 10^ years) is much shorter than the time 
required fo r these electron s to diffuse across the cool- 
ing region (jBrunettil l2003f ) . Two physical mechanisms 
have been identified as possibly responsible for the ra- 
dio emission: reacceleration of pre-existing, low-energy 
electrons in the intracluster medium (ICM) by turbu- 
lence in the core region (Gitti et al. 2002!, l2004f ). and 
the generation of secondary particles via inelastic col- 
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lisions betwee n relativistic cosmic-ray p r otons and ther- 
mal p rotons CPfro mmer fc EnfilinI 120041 : iKeshet fc Loebl 
[2012 [Keshet 2010). 

In the reacceleration model, the seed electrons may 
be provided by buoyant bubbles infiated by the central 
AGN and disrupted by gas motions in the core. In the 
absence of a reacceleration mechanism, electrons in such 
disrupted bubbles cool rapidly and emit at radio fre- 
quencies well below those currently observable. How- 
ever, a key question is the origin of the tur bulence re- 
spons ible for reaccelerating the electrons. iGitti et aTl 
(|2002f ) originally proposed that the cooling fiow of gas 
inward in the core may generate turbulence. However, 
recent X-ray observatio ns indicate that, while "cooling 
fiows" as envisioned in iFabian fc NulsenI ([^94) proba- 
bly do not materialize, even relatively relaxed clusters 
have large-scale gas motions in their cores. The obser- 
vational signature of these gas motions are the spiral- 
shaped "cold fronts" see n in the majority of cool-cor e 
clusters (for a review see iMarkevitch fc Vikhlininll2007f ). 
These cold fronts are believed to be produced by the cold 
gas o f the core "sl oshing" in the cluster's deep potential 
well. iFujita et al.l (2004) showed that sloshing motions 
can produce significant turbulence in the cluster core. 
iMazzotta fc Giacintuccil ()2008f ) discovered spatial corre- 
lations between radio mini-halo emission and cold fronts 
in the X-ray images of two clusters-the minihalos ap- 
parently contained within the region delineated by the 
cold fronts. A similar correlation is seen in the Perseus 
Cluster (Markevitch fc Churazov 2012, in preparation). 
These authors suggested the correlation arises from tur- 
bulence generated by the sloshing motions, a hypothesis 
that we test in this work. Figure [1] shows a particularly 



striking example of this in the cluster RX J 1720. 1+26, 
clearly showing that the radio emission of the minihalo 
is bounded by the cold fronts as seen in X-rays. 

In the last few years, a suite of hydrodynamic and 
magnetohydrodynamic simulations of idealized cluster 
mergers have been carried out to test the sloshing sce- 
nario for the origin of the c old fronts in cool cores 
(|Ascasibar fc Markevitchll2QQ6L hereafter AM06), deter- 
mine the physical effects of sl oshing; on the cluste r ther- 
mal gas and magnetic field (I ZuHone et al.l I2Q1QL I2Q11L 
hereafter ZMJIO and ZMLll), and make comparisons to 
observations of specifi c clusters with sloshing cold fronts 
(jRoediger et al.|[2Qllf ). For the first time, we use a high- 
resolution numerical simulation to examine the connec- 
tion between gas sloshing and radio minihalos in clus- 
ters. Our simulation models the cool, magnetized core 
gas sloshing in the gravitational potential well of a mas- 
sive, initially relaxed galaxy cluster, similar to the previ- 
ous works mentioned above. Using a filtering procedure 
similar to that emp loyed in studies of t urbulence in a cos- 
mo logical context (iDolag et al.l I2QQ5I : IVazza et al.l I2QQ6L 
I2QQ9I . I2Q1QL I2Q11L '2012), we have determined that the 
sloshing motions generate turbulence in the cluster cool 
core that is spatially contained within the volume of the 
cold fronts. 

To test the predictions of the reacceleration model for 
minihalos, we have included in the simulation a popula- 
tion of passive tracer particles to provide a Lagrangian 
description of the sloshing gas. These particles provide 
us with the ability to follow relativistic electrons associ- 
ated with a parcel of magnetized fiuid, integrating the 
change in their energies determined by the fiuid prop- 
erties (in particular, the local turbulent velocity) along 
each tracer particle's trajectory. The connection between 
turbulent velocity and the acceleration of electrons is 
treated in a "subgrid '' way- we use the prescription from 
iBrunetti fc LazarianI ()2Q07[ ). which is detailed in Section 
13.21 Given the spatial and spectral distribution of rela- 
tivistic electrons along with the magnetic field of the clus- 
ter, we produce maps of simulated synchrotron emission 
that may be compared to observations of actual mini- 
halos. Our simulation can also be exploited for evalu- 
ating the relative importance of particle reacceleration 
and cosmic ray-thermal proton collisions for the origin 
of minihalos, which will be discussed in a future paper 
(ZuHone et al 2012, in preparation). 

This paper is organized as follows: in Section 2, we 
describe the MHD simulation and the treatment of tur- 
bulence. In Section 3, we describe the treatment of rel- 
ativistic particles. In Section 4, we describe the char- 
acteristics of the turbulence that develops in the cluster 
core, and based on them, estimate the reacceleration co- 
efficients. In Section 5, we detail the evolution of the 
relativistic electrons and the resulting synchrotron and 
inverse- Compton emission. Finally, in Sections 6 and 7, 
we discuss and summarize our results and future devel- 
opments of this work. We have assumed a ffat ACDM 
cosmology with h = 0.7, ftrn = 0.3, and ftb = 0.02/i~^. 



METHODS: TREATMENT OF THE 
MAGNETIZED CLUSTER GAS 



In our simulation we solve the ideal MHD equations. 
Written in conservation form in Gaussian units, they 
are: 
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2.1. MHD Simulations 



where pth is the gas pressure, and e is the gas internal 
energy per unit volume. We assume an ideal equation of 
state with 7 = 5/3. 

We performed our simulation using FLASH 3, a 
parallel hydrodynamics/A^-body astrophysical simula- 
tion code developed at the Center for Astrophysical 
Ther monuclear Fla s hes at the U niversity of Chicago 
(Frvx ehet all 120001 : iDubev et all [2009). FLASH uses 
adaptive mesh refinement (AMR), a technique that 
places higher resolution elements of the grid only where 
they are needed. We are interested in capturing sharp 
ICM features like shocks and cold fronts accurately, as 
well as resolving the inner cores of the cluster dark matter 
halos. It is particularly important to be able to resolve 
the grid adequately in these regions. AMR allows us to 
do so without needing to have the whole grid at the same 
resolution. 

FLASH 3 solves the equations of magnetohydrody- 
namics using a directionally unsplit staggered mesh 
algorithm (USM; iLee fc Dean J 120091 ). The USM al- 
gorithm used in FLASH 3 is based on a finite- 
volume, high-order Godunov scheme combined with 
a constrained transport method (CT), which guar- 
antees that the evolved mag netic field satisfie s the 
divergence- free condition (Eva ns fc HawlevI 119881 ). In 
our simulations, the order of the USM algorithm cor- 
responds to the Piec e wise-P arabolic Method (PPM) of 
iColella fc WoodwardI (|1984l ). which is ideally suited for 
capturing shocks and contact discontinuties (such as the 
cold fronts that appear in our simulations). 

The gravitational potential on the grid is set up as the 
sum of two "rigid bodies" corresponding to the contri- 
butions to the potential from both clusters. This ap- 
proach to the modeling the potential is used for sim- 
plicity and speed over solving the Poisson equation for 
the matter distribution, and is an adequate approxima- 
tion for our purposes. The de tails of this setup may be 
found in IZuHone et al.l (|2011l ) . A comparison between 
the rigid-potential and fully-modeled self-gravitating se- 
tups for the dark matter components of merging clusters 
(and a justification of using the former method for this 
work) may be found in 7. 

We refer the reader to IZuHone et al.l ()2011l ) for the de- 
tails of our initial setup of the ICM model and the mag- 





Fig. 1. — (a) Chandra X-ray image of RXJ 1720.1+26, one of the clusters exhibiting a sloshing core and a minihalo (Mazzotta &; Giacintucci 
2008). (b) Same image with the 610 MHz radio contours overlaid (from Giacintucci et al. 2012, in preparation). Contours start from +3cr 
and are spaced by factor 2, where la = 30/iJy/beam and the beam size is 8^' x 6" . The white cross marks an unrelated head-tail radio 
source. The minihalo radio emission is confined within and traces the cold fronts that are visible as brightness edges in the X-ray image 
NW and SE of the center. 



netic field. Here we sketch only briefly the details of the 
simulation performed in this work. Our initial conditions 
consist of a massive (M -- IO^^Mq.T -- 10 keV) cool- 
core galaxy cluster, and a smaller (M ~ 2 x lO^^M©) 
gasless subcluster, separated at an initial distance d = 
3 Mpc, an initial impact parameter b = 500 kpc, and 
placed on a bound mutual orbit (the core passage occurs 
at t ~ 1.8 Gyr in our adopted time frame). The size of 
the computational domain is L = 2 Mpc, with a flnest 
cell size of Ax = 1 kpc. The maximum resolution covers 
a spherical region of r ~ 300 kpc, centered on the DM 
peak of the main cluster, which encompasses all of the 
phenomena of interest in this study. 

The ICM of the main cluster is magnetized, which 
is set up initially with a tangled magnetic field con- 
figuration with an average magnetic pressure propor- 
tional to the gas pressure, with /3pi = p/ps = 100. 
This choice is consistent with constraints put on the 
radial profiles of the magnet ic field by Faraday rota- 
tion n ieasurements JB onafed e et al. 1 120 10) and simula- 
tions CDolag et al.l lT999i : iDubois fc Tevssieni2008i ). This 
initial magnetic field configuration is not relaxed, and the 
differences in magnetic pressure and tension will gener- 
ate spurious turbulent velocities |(5v| ^ 100 km s~^ in the 
gas. To examine exclusively the turbulence generated by 
the sloshing motions, we modeled our main cluster in 
isolation for t ~ 5 Gyr, increasing the numerical vis- 
cosity during this period to damp out the motions until 
they are relatively small (|Jv| ^ 10 km s~^) and the mag- 
netic field has reached a relaxed configuration. The gas 
density and temperature profiles of the main cluster are 
largely unaffected, with the largest deviations from the 
initial profiles of ~ 10% in the central few zones, imma- 
terial for our purposes. Figure 1 of ZuHone ^et al.l (^011) 
shows the initial and final profiles of the main cluster. 

Finally, our MHD simulation includes passive tracer 
particles fiowing with the ICM. These particles will be 



used to model the effects of reacceleration and energy 
losses on relativistic electrons traveling along these tra- 
jectories fsee l3.2p . Tracer particles are stored at intervals 
of 10 Myr and from these snapshots the individual parti- 
cle trajectories are extracted. The entire simulation con- 
tains approximately 10 million passive tracer particles 
initially distributed with their number density propor- 
tional to the local gas density. 

For the details of the process of gas sloshing in galaxy 
cluster cores, e.g. how it is initiated, the formation of 
sloshing cold fronts, and its effects on the magnetic field, 
we refer the re ader to past simul ation works on slosh- 
ing (lAscasibar fc Markevitch 2006|; iRoediger et~alll2011l : 
IZuHone et al.l 12010. . .2011.) . We assume these results in 
the following discussions. 



2.2. Separating Turbulence from Bulk Motions 

We are assuming in this work that relativistic electrons 
are reaccelerated by turbulence. Therefore, in order to 
determine the acceleration rate in our simulations, we 
must determine the magnitude and location of turbu- 
lence. Since not all motions in the ICM are turbulent 
(and not all turbulent motions efficiently accelerate par- 
ticles), a distinction must be made between "bulk" or 
"laminar" flows and turbulent motions. Substructures 
and mergers will drive motions in the gas on the scale of 
the merging subclusters that are not turbulent in nature. 
This applies to our isolated binary cluster merger simu- 
lations, where the most signiflcant non- turbulent motion 
is the sloshing of the cool core, with sharp cold fronts 
and tangential velocity discontinuities. Its contribution 
to the overall velocity power spectrum needs to be prop- 
erly flltered out to examine the turbulent motions alone. 
This typically involves some process of spatial averaging, 
where the bulk motions are flltered out and only the tur- 




Fig. 2. — Schematic diagram to illustrate the method for inter- 
polating mean velocities to a point in space. Boxes represent the 
volumes over which the mean velocity is determined. The velocity 
is interpolated to any point using a triangle-shaped cloud filter. 
The shading of each box represents the approximate contribution 
of each box to the local mean. 

bulent (small-scale, random) component remains Q How- 
ever, the choice of a region to average over is often a judg- 
ment call. The simplest approaches typi cally take the 
average velocity of the cluster as a whole (Sunvaev et al.l 
l2QQ3l : lHallman fc Jelte ma 201 I j), or compute the average 
in concentric ra dial shells (^lapichino fc Niemeverl [2QQ8I : 
iLau et aTll2QQ9f ), none of which is sufficient for us in the 
presence of sloshing motions. 

One successful choice of averaging, known as the "local 
mean velocity field" technique, has been applied in both 
SPH and A MR hydrody namics simulations of galaxy 
cluste rs fe g.JPolag et al] l2QQ5l: IVazza et al.ll2QQ6L I2QQ9L 
I2Q1QI l2Qlin2Q12f ), with success in separating the turbu- 
lent and bulk velocity fields in clusters where both kinds 
of motions are significant. This technique involves divid- 
ing the computational domain into "boxes" over which 
the average velocity is obtained, and the mean velocity 
local to a point in space is interpolated from these boxes 
and subtracted from the total velocity field to yield the 
turbulent component of the velocity field. 

Following these previous works, we use the following 
simple prescription to estimate the turbulence in our sim- 
ulations. We begin by setting up a coarse 3D uniform 
grid of cubical boxes of width £ that spans the entire 
simulation domain. Within each box, a mean velocity 
is computed from the cells on the FLASH adaptive grid 
contained in that box. At any point in the simulation, 
a local mean velocity field v(x) may be constructed by 
interpolation from these boxes vi a a "Triangle- Shaped 
Clou d" (TSC) window function (jHocknev fc Eastwood! 
I1988L a schematic illustration of this procedure is given 

'^ A more sophisticated approach involves the implementation of 
a subgrid model for the unres ol ved turbulent motions , such as in 
[Scannapieco fc BriiggenI (|2008l ) ; llapichino et al.l (|201in 



in Figure [2j). From this local mean field, we may compute 
a local "filtered" turbulent field Jv = v — v. 

The resulting turbulent velocity field is a conserva- 
tive estimate of the field that may arise in a real cluster 
merger such as the one we are simulating. The TSC inter- 
polation will filter out turbulent modes with scales larger 
than the box size £. Additionally, the dissipation scale 
where the energy in the turbulent cascade is dissipated 
into heat is set not by the actual viscosity of the ICM 
but by the resolution of our simulation, which results in 
a dissipation scale of several times the cell size. Thus, the 
power of turbulence that we measure in our simulations 
is a lower bound to that which may be expected in a real 
cluster merger. We will address this point in detail in 
Section 14.11 We will also discuss a separation of turbu- 
lent motions into compressive and solenoidal components 
in Section [ 



3. METHODS: TREATMENT OF RELATIVISTIC 
ELECTRONS 

The reacceleration model for radio halos assumes a 
population of relativistic electrons already present in the 
ICM, which are reaccelerated by turbulence, e.g. that 
generated by sloshing motions. In our simulation model, 
we assume that relativistic electrons are advected along 
with the fiuid, so we may use our tracer particle tra- 
jectories to determine the effect of the properties of the 
fluid along each trajectory on the energies of the rel- 
ativistic electrons. We do this by performing a set of 
post-processing "simulations" on the tracer particle tra- 
jectories. There is a suflicient number of tracer parti- 
cles to sample the cool core very well. A similar post- 
processing approach was used to study simulated gas 
mixing, metal distribution, and iron line profiles in a 
set of cluste r s dra wn from a cosmological simulation in 
IVazza et all (J2Q1Q[ ). 

3.1. Model for Relativistic Electron Spectra 

In order to compute the synchrotron emission for each 
tracer particle as a function of time, it is necessary to 
model the energy distribution of relativistic electrons 
A^(7) that may be initialized and evolved for each tracer 
particle along its trajectory. A rigorous approach to this 
problem would be to integrate a Fokker-Planck equa- 
tion for N{j) along each trajectory, with N{j) sam- 
pled discretely in bins of electron energy. This ap- 
proach would provide a complete description of the rel- 
ativistic particle diffusion in momentum space as well 
as systematic energy gain and losses. While feasible 
in the near future, this approach is beyond the cur- 
rent state of software development for our code, and we 
chose to adopt a simpler method. Instead, the distri- 
bution of relativistic electrons A^(7) is discretely sam- 
pled in a "Monte Carlo" fashion with a large number 
of "relativistic particles" for each tracer particle. The 
distinction between these two methods is similar to the 
distinction between Eulerian grid-based methods versus 
Lagrangian particle-based methods (smooth-particle hy- 
drodynamics) for solving the equations of hydrodynam- 
ics. In the latter case, the individual gas particles do not 
represent any real concentrations of gas per se^ but are 
simply discrete Lagrangian samples of the distribution 
function of the gas; similarly, the relativistic particles 



are not electrons but merely sample the underlying dis- 
tribution. 

Each relativist ic particle has a Lorentz factor 7^ j that 
represents the average relativist ic energy jijrrieC^ for the 
sample i for the tracer particle j. The energy disribution 
function Nj (7) for each tracer particle j may be written 
as: 



iV,(7)=if,-V.5(7-7i,,), 



(7) 



where the sum is taken over the number of samples for 
the tracer particle and Kj is a normalization constant. 
To get a proper scaling for the total energy in relativistic 
electrons, as our default case we assume the initial energy 
density in relativistic electrons of each tracer particle is 
proportional to the tracer particle's associated thermal 
energy (that is, assuming the uniform initial ratio of rel- 
ativistic electron to thermal pressures): 
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where £j is the internal energy per unit mass of the ther- 
mal gas and each tracer particle has a gas mass rrij asso- 
ciated with it, such that the total gas mass of particles 
equals the total mass in gas on the grid within the region 
under consideration. This condition sets the normaliza- 
tion constant Kj. The results in this work only depend 
on 77 as a normalization constant, so we will express elec- 
tron energies and radiative quantities in terms of 77, but 
we assume a default value of 77 = 10~^ to compare the 
results of our numerical simulations to observations. 

We generate a physically plausible initial spectrum of 
relativistic electrons by taking a power-law spectrum 
with index s = 2.5 (defined as A^(7) ex 7"*) and pas- 
sively evolving it for 1 Gyr under the constant conditions 
of 71th = 0.01 cm~^ and B = 5 /iG, representative values 
for the cluster core regions in this model. The result is 
a spectrum that has significantly cooled, with essentially 
no electrons at energies higher than 7 ^ 500. This initial 
distribution of electrons has no significant emission at 
observable radio wavelengths. We assign this spectrum 
to each tracer particle at t = 2.55 Gyr. At this epoch, 
the characteristic size of the sloshing region is approxi- 
mately 41 = 120 kpc, and by beginning the relativistic 
electron evolution at this point, we hope to minimize 
contamination to our derived turbulent spectrum from 
bulk motions. This is also roughly the epoch at which 
turbulence begins to be significant (Section 14. ip . 

As a starting point, we also assume that the number 
density of relativistic electrons is proportional to the den- 
sity of thermal electrons (however, in Section [5^ we will 
discuss alternative initial spatial distributions of parti- 
cles). Since our tracer particles are advected with the 
thermal gas and represent parcels of constant gas mass, 
this simply implies that the number of relativistic elec- 
trons per tracer particle remains constant. For the re- 
sults presented here, we generated 10^ relativistic parti- 
cles for each tracer particle. 

We select only certain tracer particles for consideration 
instead of using all of them contained within the simu- 
lation volume. Our default case begins with the tracer 
particles filling a large region of radius R = 300 kpc in 



the initial undisturbed cluster, covering the region that 
will eventually contain sloshing and cold fronts. 

3.2. Two Mechanisms of Reacceleration by MHD 
turbulence 

The critical component of the model for radio mini- 
halos is the reacceleration of relativistic particles by 
MHD turbulence. Relativistic particles interacting with 
MHD turbulence are subject to several mechanisms that 
lead to their acceleration, both through resonant and 
non-resonant coupling with turbulent waves. In the 
case of compressible MHD turbulence generated at large 
scales in the ICM, B runetti fc Lazarian (2007 ) have 
shown that reacceleration is accomplished primarily by 
the damping of fast magnetosonic waves thro ugh the 
Transit-Time Damping (TTD) resonance (e.g., iMelrosd 
119681) . the condition for which is 



k\ivi\ =0, 



(9) 



where u is the frequency of the magnetosonic wave, 
and /c|| and v\\ are the parallel (projected along B) 
wavenumber and particle speed, respectively. To com- 
pute the momentum-diffusion coefficient of relativistic 
electro ns due to TTD resonance, we use Equation 40 
from .Brunetti fc Lazarian (2007.) (assuming the mode 
phase velocity (Iph) ^ c^), where the spectrum of com- 
pressive modes is calculated self-consistently by taking 
into account the effect of dampings due to TTD with 
thermal and relativistic particles: 
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In the integral over /c, /cmin specifies the wavenumber 
corresponding to the scale of the largest turbulent ed- 
dies and /ccut specifies the wavenumber corresponding to 
the minimum scale of turbulence where the timescale of 
dampings equals that of the turbulent cascade, p is the 
electron momentum, B^ is the local magnetic field, H 
is the Heaviside step function, and the spectrum of the 
magnetic field fiuctuations associated with the compres- 
sive modes >Vb(/^) is 
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(11) 



where the quantity in brackets (...) is averaged with re- 
spect to Q and is of order unity. Noting that /3pi = 
Pth/pB^ Equation [To] may be simplified to 
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where X is the integral over and is of order a few. The 
integral over k may be rearranged to give 

dkW{k)k = (k) [ '"' dkW{k) ^ {k)pvt^R'', (13) 

where {k) is the average wavenumber of the turbulent 
cascade as measured in the simulation, Vt is the turbu- 
lent velocity of the tracer particle, and R^ is the frac- 
tion of the turbulent energy in the form of compressible 
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motions. The last equality follows from the definition 
of the power spectrum. The scaling of fast modes with 
wavenumber k is of Kraichnan form, with W(/c) ex /c~^/^ 
(JCho & Lazarian 2003), as opposed to the purely hydro- 
dynamic Kolmogorov scaling of k~^^^. We will assume 
this form for the power spectrum of the turbulent cas- 
cade throughout this paper. The mean wavenumber (k) 
is then expressed in terms of /cmin and kcut as: 



(k) 



"^cut 



W{k)kdk I "' k-^^^kdk 



I^cut 



W{k)dk 



"^cut 



V "^min'^c 



fc-3/2dfc 



(14) 
In cgs units, the resulting momentum-diffusion coeffi- 
cient for the relativistic electrons via TTD is given by: 



D, 



pp,TTD 



4x 10--''(fc)vt i^V 



(15) 



Fast particles diffusing in compressive MHD turbulence 
also experience statistical compression and expansion 
that leads to additional stochastic non-resonant accelera- 
tion. In the f ast diffusion limit, the mon ientum- diffusion 
coefficient is (jBrunetti fc LazarianI 120071 Equation 55): 
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where the spatial di ffusion coefficient D under our physi- 
cal conditions is is (jBrunetti fc Lazarianll2007l Equation 

56): 
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The two mechanisms, TTD resonance and non-resonant 
compression, are essentially driven by the same turbu- 
lent modes and involve independent particle-mode cou- 
plings. Consequently, as a first approximation, the ac- 
celeration process may be tho ught of as the combi- 
nation of the two mechan isms (jCho fc LazarianI 120061 : 
iBrunetti fc La zarian' '2007), with a total momentum- 
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3.3. Gains and Losses on the Relativistic Particles 
Each relativistic particle has a reacceleration efficiency 
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SO that the rate of momentum increase for each relativis 
tic particle is 
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We will detail the assumptions that go into the compu- 
tation of the reacceleration coefficient in our simulations 
in Section 1421 ^^ particular, the evaluation of (/c), i?^, 
and Vt. 

The rate of energy losses due to synchrotron and 
Inverse- Compton scattering off of cosmic microwave 
background (CMB) photons for individual electrons is 
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Since our relativistic particles are not individual elec- 
trons but samples of the electron distribution function 
(and hence represent many electrons), we assume for 
simplicity that each relativistic particle represents an 
isotropic distribution of pitch angles, with (sin^ 0) = 2/3. 
^CMB = 3.2(1 + zYiiG is the equivalent magnetic field 
strength for the CMB at present, where z is the cosmo- 
logical redshift. 

In our simulation the core passage of the disturbing 
subcluster occurs at t ~ 1.8 Gyr from the beginning of 
the simulation, and our particle trajectories begin at t 
= 2.55 Gyr (after the onset of turbulence due to slosh- 
ing, see Section [5TT]) . and we assign the redshift z = 
to the epoch t = 5 Gyr of the simulation, in order to 
reproduce some of the observed nearby clusters exhibit- 
ing cold fronts in their cores. Under these conditions, 
the redshift at each epoch is computed from the simula- 
tion time assuming a /i = 0.7, VLm = 0.3, and Qa = 0.7 
ACDM cosmology. 

The Coulomb losses are given by (in cgs units): 
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where nth is the number density of thermal particles. 

3.4. Solving for the Evolution of the Relativistic 
Particles 

If spatial diffusion is not important, formally the time 
evolution of the relativistic electron momentum distribu- 
tion N(p,t) is a solution to the Fokker-Planck equation 
(jBrunetti fc Lazarianll2007[ ): 
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Solving this equation numerically can be expensive, par- 
ticularly for the case of many individual tracer particle 
trajectories as in our simulation. However, as previously 
described, we have chosen to evolve relativistic "sam- 
ple" particles instead of N{p^ t) explicitly, which can be 
thought of as the probability density for the random 
variable P^, which corresponds to the momenta of the 
sample particles. For a given Fokker-Planck equation 
and distribution function N{p^ t) there is a correspond- 
ing stochastic differential equation (SDE) for the evo- 
lution of Pt for an ensemble of sample particles. By 
following the momentum trajectories of many sample 
particles, we may reliably approximate the behavior of 
N{p^t) and the observable quantities that depend on 
it, such as the resultant synchrotron and IC emission. 
SDEs have been used extensively in other astrophysi- 
cal contexts, in particular, for the integration of cosmic 
ray trajectories in the heliospheric and galactic magnetic 



fields (IZhandllQQalFlQm fc PQg:QrelQvl2QQ9l : [Pei et all 



I2Q1QI : iStrauss et all 1201 ll : iKopp et alJl2Q12f ). The main 
differences between our approach and many of these 
works is that a) we are only integrating the momentum 
of each particle as the relativistic particles are assumed 
to follow the tracer particle trajectories in space, and b) 
we are integrating the equations forward in time instead 
of backwards to the original source of particles. 

The SDE that corresponds to the above Fokker-Planck 
equation is given (in the Ito formulation) by: 



dPt = a{p, t)dt + h{p, t)dWt 



where the "drift" term is 
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and the "stochastic" term is 



b{p, t)dWt = y^2D^^dWt - ^2DppdtM{^, 1) (25) 

where dWt is a standard Wiener (or " Brownian motion" ) 
process, and A/'(0, 1) is a normal distribution with zero 
mean and unit variance (the "~" symbol here indicates 
"is distributed as"). The effects of each of these terms 
on the relativistic electron spectrum will be shown in the 
Appendix. 

We have integrated the drift term of this equa- 
tion using a 4th-order Runge-Kutta method. To inte- 
grate the stochastic term , we use the Milstein method 
(jKloeden fc Platen! 120111 ). This results in the following 
discretization for Equation [25l 

h{p, t)dWt « h(p, i)AVr„ + h,(p, i)^^M[(AM^„)2 _ At] 

(26) 
where 

AWn - V^Ar(0, 1) (27) 

This equation is integrated for each relativistic par- 
ticle along each tracer particle trajectory with a 
variable timestep for each tracer particle /Sij = 
min{0.1(pi,j/pi,j),0.1(pf .^/I^pp^ij)} to ensure stability. 
To determine the fluid quantities (p, T, 5v^ B) at any 
point on the trajectory, they are linearly interpolated 
between the saved instances of each particle along that 
trajectory, which have a time resolution of 10 Myr. In the 
Appendix we provide tests of our method when compared 
to analytical solutions and the results of a Fokker-Planck 
calculation. A detailed description of the method and the 
code will be the subject of a future paper (ZuHone et al 
2012, in preparation). 

4. RESULTS: TURBULENCE 

4.1. Characteristics of Turbulence Generated in the 
Sloshing Region 

In the sloshing/reacceleration hypothesis for radio 
mini-halos that we are testing in this work, the radio 
emission coincident with the envelope of the cold fronts 
as seen in X-rays is due to the turbulence that is associ- 
ated with the sloshing motions. Figure [3] shows the spiral 
shape in temperature that results from the sloshing mo- 
tions for a few different epochs of the simulation. We 
now determine the location and spectrum of turbulence 
that results from the encounter with the subcluster. We 
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Fig. 3. — Projected temperature maps in the 2;-projection for 
several different epochs of the simulation. The colorbar is temper- 
ature in keV. Each panel is 750 kpc on a side. Tick marks indicate 
a distance of 100 kpc. Core passage occurs at t approxl.8 Gyr. 



note at the outset that a very recent work (Vazza e t al.l 
2012), using an improved filtering technique, produced 
turbulent velocity maps and power spectra for gas slosh- 
ing in the core of a galaxy cluster and achieved results 
in general agreement with those we present here. 

Separating out the turbulent energy component from 
the bulk component in a simulation of gas sloshing is a 
challenge, mainly due to the presence of the sloshing mo- 
tions themselves. Strong velocity shears and discontinu- 
ities exist at the cold front surfaces, and tracer particles 
situated near these surfaces will have a mean velocity 
interpolated to them that is a weighted sum of the ve- 
locities across the front, resulting in some of the bulk 
component of the particle's velocity remaining after the 
spatial filtering procedure. Therefore, tracer particles 
that travel near these surfaces may have their turbulent 
kinetic energy overestimated, which we should find a way 
to avoid. 

In our merger setup, where the subcluster's trajectory 
is in the x — y plane, most of the kinetic energy associated 
with the sloshing motions is in the Vx and Vy components 
of the velocity field. It is in these directions in which it 
will be most difficult to filter out the contribution to the 
total velocity at a particle position from the sloshing mo- 
tions. At the same time, we expect that the vast majority 
of the kinetic energy associated with the Vz velocity com- 
ponent is due to the turbulent cascade generated by the 
sloshing motions, since there is only a slow outward ex- 
pansion of the cold fronts in the z-direction. We will use 
this simple observation to separate the turbulent compo- 
nent in our simulation. 

To get a sense of the location of the tubulence in 
our simulations, we create projected maps of the mass- 
weighted turbulent velocity: 
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where the integrals are taken along the line of sight. The 
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Fig. 4. — Projected (mass- weighted) turbulent velocity map using all velocity components in km s"-*^ at the epoch t = 3.35 Gyr for 
varying filtering box size i. Left panels: Projection of all three velocity components. Right panels: Projection of the Vz component only. 
In each panel a sample filtering box is shown. Each panel is 750 kpc on a side. Tick marks indicate 100 kpc distances. 



left panel of Figure |4] shows the projected mass- weighted 
total velocity at the epoch t = 3.45 Gyr, ~1.7 Gyr after 
the subcluster core passage, for a few different choices 
of the filtering box size i. There is a clear distinction 
between the core region bounded by the sloshing cold 
fronts and the region outside; inside this region, (5vrms ^ 
50 — 200 km s~^, while outside the cold fronts, Svrms ^ 
10 — 20 km s~^. However, the highest values of projected 
J'Urms appear to be near the cold front surfaces, which 
is to be expected from the preceding discussion. The 
speed of these motions increases with increasing box size 
^, as the larger box size permits larger-scale flows to be 
interpreted as "turbulent" by the filtering scheme. 

The right panel of Figure |1] shows the projected mass- 
weighted z-component of velocity at the same epoch 
t = 3.45 Gyr, for the same choices of the filtering box 
size i. In this case, there is still a clear distinction be- 
tween the sloshing core region and the outside region of 
the simulation domain. However, the Vz field does not 
in general have anomalously large values along the cold 
front surfaces. This indicates that (given our filtering 
procedure) the Vz component of the filtered velocity is 
a better indicator of the turbulent pattern in our sim- 
ulated ICM. As expected, the turbulent energy within 
the sloshing region is smaller than in the case where the 
Vx and Vy components of velocity were included. In our 
simulations, we will use dvz to estimate the turbulent 
velocity, scaling Svz by an appropriate factor to reflect 
the fact that the total kinetic energy in turbulence is 
larger than just the contribution from this component 
alone. Note that because it is the sloshing and shear in 
the X — y plane that drives turbulence, the true Vx and 
Vy turbulent components will be relatively higher than 
they would be in isotropic turbulence, and this should 
be taken into account. We also expect, on the basis of 
Figure HI that the overall spatial pattern of turbulence 
in all three directions will be similar, so by taking the 
Vz component only, there is no significant loss of spatial 
information. 

It is instructive to compare the power spectra of the 
filtered and unfiltered velocity field, as well as the spectra 
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Fig. 5. — The total velocity power spectrum for differing filtering 
box sizes £. The epoch is t = 3.45 Gyr. 



of the velocity field for different filtering lengths. For 
this, we wish to examine the region that is dominated 
by the sloshing motions, so the power spectrum is taken 
over a box size of L = 250 kpc on a side, centered on 
the potential minimum of the cool-core cluster. The 1-D 
power spectrum P{k) of the velocity field is defined such 
that 



(«2) = :^||v(x)pd3^ 



27r2y 



P{k)k'^dk 



(29) 



which follows from Parseval's theorem. 

Figure [5] shows the total velocity power spectra for a 
few different values of the filtering box scale l^ in com- 
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Fig. 6. — The velocity power spectrum for the different compo- 
nents of velocity for the filtered (dashed lines) and unfiltered (solid 
lines) components of velocity. For the filtered components, the 
filtering box size is ^ = 30 kpc. The epoch is t = 3.45 Gyr. 



parison to the unfiltered power spectrum at the epoch 
t = 3.45 Gyr of the simulation. The unfiltered power 
spectrum maintains a rough power-law shape from low to 
high wave number. This indicates that the sloshing mo- 
tions contribute power at nearly all scales (as expected, 
given the geometry of the fronts with sharp velocity dis- 
continuities). In contrast, the filtered power spectra are 
characterized by the following features: 1) an increase 
at low wavenumbers, followed by 2) a roughly power- 
law scaling at intermediate wavenumbers, and finally 3) 
a drop-off at high wavenumbers. The first part of the 
spectrum can be explained by the smoothing out of ve- 
locity structures with scales larger than the filtering box 
size i due to our filtering procedure. The second part of 
the power spectrum corresponds to an "inertial range" 
with a power-law scaling P{k) ex /c~". The third part of 
the power spectrum, a steep decrease in power beginning 
at wavenumbers corresponding to linear scales ^ 8 Ax, is 
due to the effect of the numerical dissipation associated 
with the PPM hydro method used by our si mulations (as 
noted by [Por ter & Woodward 1994; Vazz a et al.l I2QQ9I : 
iKitsionas eFa l. 2009; Vazza et al.. .201 1)1^ The unfortu- 
nate side effect of the filtering box is that it filters out 
contributions at low wavenumbers that may be part of 
the turbulent cascade. As i is decreased, the normaliza- 
tion of the power spectrum also decreases. However, we 
find that that for filtering box sizes £ < 30, the power 
spectrum converges at intermediate to high wavenum- 
bers, indicating that at these scales, we have largely fil- 
tered out the sloshing motions. 

Figure [6] shows the power spectrum in the three differ- 
ent velocity components, for the filtered (with a filtering 
box size of ^ = 30 kpc) and unfiltered velocity fields, at 
the same epoch t = 3.45 Gyr. Though the curves are 
roughly the same shape, the power spectra of the x and 

^ Since this is a dissipative effect on the small-scale turbulent 
motions and not on the sharp velocity discontinuities of the cold 
fronts, we do not see it in the unfiltered spectrum. 



y components of velocity are each a factor of ~2-3 larger 
than the power spectrum of the 2:-component of velocity 
for all /c, a reflection of the fact that sloshing motions 
which are responsible for the generation of turbulence 
are mostly in the x — y plane. This indicates that the to- 
tal kinetic energy is roughly ~5-6 times higher than the 
kinetic energy in the 5vz component of velocity alone 
(compared to the factor of 3 for isotropic motions). 

The left panel of Figure [7] shows the power spectrum 
at varying epochs, with a filtering box size ^ = 30 kpc. 
At the beginning of the simulation, prior to the subclus- 
ter infall, there is a small degree of turbulence, due to 
the spurious velocities generated during the relaxation 
of the initial condition (Section 12. ip . Once the slosh- 
ing begins, after the core passage that occurs around 
t ^1.8 Gyr, the power spectrum normalization increases 
until approximately t = 3 Gyr, after which the normal- 
ization is roughly constant. This is due to the fact that 
at these late epochs, the volume of our 250 kpc box no 
longer contains all of the sloshing motions, which con- 
tinue to develop turbulence outside this region. Addi- 
tionally, we note that the shape of the power spectrum 
curves are very similar between all epochs. The right 
panel of Figure [71 shows the power spectrum at the same 
epochs, multiplied by /c^/^. Within the inertial range of 
wavenumbers, this quantity at late epochs is fairly flat, 
indicating rough agreement with our assumption that the 
power spectrum is of Kraichnan {P{k) ex /c~^/^) form at 
these wavenumbers (Equations [TQlfT9l) Fl 

4.2. Estimating Reacceleration Coefficients 

We may use the information regarding the spectrum 
and spatial distribution of turbulence from our simula- 
tion to estimate the reacceleration coefficient at the posi- 
tions of the tracer particles (Equations [TQlfT9l) . As men- 
tioned above, our procedure for estimating turbulence 
should account for a) our use of only the Vz component 
of velocity to minimize the influence of non-turbulent mo- 
tions in the velocity field, b) the fact that only compres- 
sive MHD turbulence reaccelerates relativistic particles 
efficiently, and c) the turbulent cascade being artificially 
cut off at large scales by the size of the filtering box and 
at small scales by the resolution scale of the simulation. 
The first and third of these effects result in an underes- 
timate in the total turbulent kinetic energy available for 
reacceleration, and the second results in an overestimate. 
The following represents an attempt to account for these 
effects in an averaged way. 

The first consideration is the underestimate in the ki- 
netic energy due to only taking the Vz component. We 
have already determined that the power spectra of the 
three different components have roughly the same shape 
(see Figure [6]), and the total turbulent kinetic energy in 
all three components is roughly a factor of 6 larger than 
the Vz component alone, which we will adopt in the fol- 
lowing estimates Q. 

^ Adopting a Kolmogorov power spectrum (P(k) oc k~^^^) 
changes the numerical coefficients only slightly and does not have 
a substantive effect on our results. 

^ Formally, Equation [To] assumes isotropic fast modes with re- 
spect to the local magnetic field direction. Still, we adopt the 
form of Equation 40 from Brunetti & Lazarian (2007), noting that 
checking whether this condition is completely satisfied is beyond 
the aim of this paper 
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Fig. 7. — Evolution of the 2;-component velocity power spectrum vs time. The filtering box size is ^ = 30 kpc. Left panel: The power 
spectrum P{k). Right panel: The scaled power spectrum k^^'^P{k). 



The second consideration is regarding the fraction of 
turbulent energy in the form of compressive modes, R^^ 
that we use for the acceleration of particles in our par- 
ticular model. The turbulence in our case is composed 
of both solenoidal and compressive components; a num- 
ber of previous investigations have shown that for a 
purely solenoidal driving force, that the resulting tur- 
bulent cascade will be l argely incompressible . For hy- 
drodynamic turbulence, ("Bert oglio et al.|[2QQl[ ) estimate 
a ratio Eq/Es cx Re x M^, Mt and Re being the tur- 
bulent Mach number and the Reynolds number in the 
medium (though this scaling only applies if Re x M^ < 
10). Using numerical simulations of MHD turbulence, 
iCho fc LazarianI (|2QQ3f ) showed that the drain of energy 
from solenoidal to compressive motions is fairly smah 
for subsonic (and super-Alfvenic) turbulence, consistent 
with (their Equation 1): 
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where Va is the Alfven speed and {SV)a is the charac- 
teristic velocity of the Alfvenic turbulent modes. Conse- 
quently, under our subsonic, moderately super-Alfvenic 
condition^, we expect that about 10% (or less) of the 
energy of incompressible motions is transferred into the 
compressive cascade, simply meaning that if both com- 
pressible and incompressible modes exist they generate 
independent cascades. 

For subsonic turbu lence, R^ is essentially determined 
by the driving force (jFederrath et al.|[2Qlll ). In our sit- 
uation, where two dark matter cores pass by each other 
with a small impact parameter, both compressible and 
solenoidal driving is present in the region of the slosh- 
ing motions. We attempt to estimate the fraction of 
compressive turbulent power in our sloshing core as fol- 
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lows. Following oth er authors (e.g., iRvu et ahl I2QQ8I : 
iFederrath et al .1120 111 ) , we estimate the ratio of compres- 
sive to total power R^ in our simulation by applying a 
decomposition in the Fourier space of the velocity field 
into its compressive and solenoidal components. For each 
of these components, we may compute the power spec- 
trum which will yield an estimate of R"^ for a given vol- 
ume. This quantity cannot be computed at each point 
in space; therefore, exactly which volume to choose is 
an important consideration, because there is a range of 
turbulent velocities at any given epoch within the core, 
which may be associated with a range of values for R^. 
Since our filtering procedure essentially filters out con- 
tributions to the velocity field on lengths larger than the 
filtering box size £^ this sets a natural scale that should be 
sufficient to provide a " local" estimate of R^ for a range 
of turbulent velocities within the cluster core region. 

For several epochs of the simulation, we generate uni- 
formly gridded velocity data (at our finest resolution of 
Ax = 1 kpc) within a box size of L^ = (300 kpc)^ 
centered on the cluster potential minimum, and divide 
this data into 10^ smaller boxes of size i^ = (30 kpc)^ 
each, the size of our filtering boxes. For each of these 
small boxes, we take the Fourier transform of the veloc- 
ity components Vi (x) , and separate the compressive and 
solenoidal components of the transformed velocity field 
Vi{k.) by the following projection operations in /c-space: 



\k) = kikjVj{k) 
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(32) 



where the indices (i, j) indicate the different spatial com- 
ponents and the Einstein summation convention over re- 
peated indices is assumed. For each of these components 
a power spectrum may be computed in the usual manner 
f Equation [29]) . 

Figure [8] shows phase plots of the fraction oi V = l^ 
domains with a given average 3D turbulent velocity Vt 
and ratio R'^ = Pc/Ptot over several epochs of the sim- 
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Fig. 8. — Phase plots of turbulent velocity vt and the ratio R^ 
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ulation within the central L^ volume. At t = 2.55 Gyr 
(0.75 Gyr after core passage), near the beginning of the 
sloshing period, the turbulent velocities are mostly very 
low {vt^bO km s~^-most of the sloshing motions are 
bulk flows), with ratios of R^ clustering around ~0.1- 
0.3, and a tail of points extending up to R^ ~ 0.6. As 
the sloshing motions expand and drive turbulence, the 
distribution of turbulent velocities spreads out over the 
range ~0-250 km s~^, clustering mainly within the ~50- 
150 km s~^ range, until the epoch t = 3.95 Gyr, after 
which the average velocity in the majority of these re- 
gions steadily decreases to the ~50-150 km s~^ range, 
though a tail of strongly turbulent regions always re- 
mains. Throughout this entire time period, the value of 
R"^ for most of these regions stays within ~0. 1-0.3, with 
a tail of regions extending up to i?^ ~ 0.6, and a typical 
value of R'^ ~ 0.25 for cells with higher {vt ^ 100 km s~ ) 
turbulent velocities, where acceleration is important. In- 
terestingly, the region of the phase plots with high Vt 
and high R^ is mostly devoid of points, likely related to 
the fact that the numerical viscosity inherent in the sim- 



ulation preferentally damps compressible motions. As a 
consistency check, we compute the ratio Re for the entire 
region dominated by the sloshing motions for the same 
epochs, and we find that during this time Re ~ 0.25 — 0.5, 
consistent with our more spatially resolved estimates. 

We conservatively adopt the value i?^ = 0.25 as the de- 
fault value for our calculations of the reacceleration and 
momentum-diffusion coefficients. Adopting the average 
value is conservative, because what is more relevant in a 
non-linear process such as particle reacceleration is the 
tail of high values of R^ and Vt^ not their average. 

A caveat must be made regarding this procedure for 
determining R^. The Fourier transform over a finite do- 
main assumes that the field is periodic. Computing the 
Fourier transform of a velocity field on a non-periodic do- 
main is equivalent to taking the transform of a field that 
has a sharp discontinuity such as a shock or a cold front 
at the boundaries of the domain. This will add spurious, 
unphysical power into the computation in both the com- 
pressive and solenoidal components of the velocity field, 
though the amount of spurious power should be small if 
the domain over which the FFT is taken is larger than 
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the largest scale at which significant power is present in 
the velocity field. 

In order to make a rough determination of the spuri- 
ous compressive power that this effect will introduce, we 
have performed the same analysis on a velocity field that 
has no compressive power by construction. To construct 
this field, we initialize a Gaussian random velocity field in 
/c-space on a L^ = (300 kpc)^ domain using a power spec- 
trum very similar to that found in our simulation, and 
perform a divergence-cleaning operation on this field to 
remove the compressible component from the field. This 
field is Fourier-transformed to real space, and the same 
analysis is performed on this field as on the velocity field 
from our cluster simulation. Since there is no compres- 
sive power in this field, any Re > in the smaller {i = 
30 kpc) boxes must arise from the aforementioned effect 
at the edges due to the non-periodic field in the subdo- 
mains (we have verified that over the entire domain we 
find Re = to machine precision, since it is periodic). 
We find that the value of Re in these boxes ranges from 
0-0.1, with a mean value of Re ~ 0.07. This indicates 
that about this much of the compressive power that we 
estimate from our sloshing simulation is potentially spu- 
rious. 

The third consideration regards the inertial range 
of the power spectrum and its effect on the average 
wavenumber (k) as well as on the total estimated ki- 
netic energy of turbulence. The power spectrum begins 
to drop off from the inertial range at high wavenumbers 
due to the dissipation associated with the finite resolu- 
tion of the simulation. In reality, there will be a physical 
damping scale, and we should try to use that scale, and 
not the artificial numerical one, for calculating the reac- 
cleration coefficient. 

Turbulence under conditions in cluster cool cores is 
mostly collisional. Collisionless damping of the magne- 
tosonic waves with the thermal plasma becomes strong 
as soon as turbulence reaches the electron Coulomb 
mean free path -^mfp, and it is at scales similar to 
this scale that we expect t he turbulent cascade to end 
(jBrunetti fc Lazarianl '200"7) . If plasma instabilities play 
a role, they will make the ICM more collisional, and we 
may expect that the turbulen t cascade would exte nd to 
much smaller scales ( Brunetti fc Laz arianI l2011b[ ). To 
be conservative, we adopt the standard picture where 
the turbulent cascade is cut off by collisionless damp- 
ing on thermal particles. For the conditions in the 
core of our model cluster, the electron mean free path 
^mfp ^ 0.2 — 0.1 kpc. If we assume that the iner- 
tial range of the turbulent cascade extends to at least 
kcut ^ 27r/^mfp, we conservatively estimate that the ki- 
netic energy in turbulence should be increased, over what 
we obtain by simply integrating the spectrum in the sim- 
ulations, by a factor of / ~ 1.5. 

Taking these considerations into account, we make the 
following modifications to Equation [15] for the TTD co- 
efficient: 
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Fig. 9. — Relativist ic electron spectra for several different epochs 
of the simulation, beginning with the initially injected spectrum at 
t = 2.55 Gyr. 



efficient: 



l^pp,c - 1.3 X 10-^^fe. 
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(34) 

where v'^/v'^ ^ 6, /cmfp = 27r//mfp, and the approxi- 
mate correction for extending the inertial range to high 
wavenumber is / ~ 1.5. The assumptions regarding the 
extent of the power spectrum also affect the computation 
of the average wavenumber (k). Assuming (/cmin,^cut) = 
(27r/30 kpc"\27r/0.1 kpc"^), from Equation [ll we find 
(k) ^ 3.6 kpc~^. For our conditions, the coefficients 
^pp,TTD and ^pp,c are of a similar order of magnitude, 
and they are added to produce the total ^pp. The var- 
ious uncertainties associated with these corrections are 
discussed in Section [ 



(33) 
and combining Equations flBlflTl for the non-resonant co- 



5. RESULTS: THE EVOLUTION OF RELATIVISTIC 
ELECTRONS, SYNCHROTRON, AND IC 
EMISSION 

5.1. The Evolution of Electron Spectra 

Though unobservable directly in real clusters, from 
our simulations we may examine the relativistic elec- 
tron spectrum of the cluster as a function of time dur- 
ing the simulation. We construct the energy spectrum 
^(7) by binning up the relativistic particle samples into 
100 equally log-spaced bins over the range (7min, 7max) = 
(10, 10^). For simplicity, we include all of the relativis- 
tic particles in the cluster in the binning procedure (we 
will examine the spatial distribution of the relativistic 
electrons in the next section). 

Figure [9] shows the relativistic electron energy spec- 
trum N{j) for all tracer particles for several different 
epochs of the simulation, with times given from the be- 
ginning of the simulation. Only 0.2 Gyr after the injec- 
tion of relativistic particles, reacceleration has already 
generated a population of particles with energies up to 
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Fig. 10. — Projected gas temperature maps with radio contours overlaid at several epochs for the frequencies 153, 327, and 1420 MHz in the 
2;-projection. The colorbar is temperature in keV. Contours of radio emission at (153, 327, 1420) MHz begin at (1.0,0.5,0.125) x 10""^ mJy 
arcsec""^ and increase by a factor of 2. Each panel is 750 kpc on a side. Tick marks indicate a distance of 100 kpc. 



7 ^ 10^0 At later times (over the course of approxi- 
mately 2 Gyr), reacceleration maintains a population of 
relativistic electrons up to 7 ~ 2 x 10^. In our particular 
merger setup, reacceleration can no longer keep up with 
the cooling of the relativistic particles about 1.5 Gyr af- 
ter the core passage (t ~ 3.5 Gyr), and the spectrum 
begins to steepen. 

5.2. Simulated Synchrotron Radiation 

Once we have the accelerated electron spectra at a 
given epoch for each tracer particle, we can use them 
to compute the synchrotron radiation they emit. The 
synchrotron power for a single electron as a function of 
frequency is (jRybicki fc Lightmanlll979f ): 



P{r^,l) = ^^^Fi.) 



rrieC^ 



(35) 



where F{x) is the synchrotron function, and x = ly/i^c^ 
where Uc = {3/4:7T)j^eB±/rrieC is the synchrotron criti- 
cal frequency. The total synchrotron power at a given 
frequency for each tracer particle is then the sum of syn- 
chrotron powers for the electron samples associated with 
the tracer particle scaled by the normalization constant: 



li: 



(36) 



i 

^ Since there is such a large change in the electron spectrum 
within the first 0.2 Gyr, we do not expect much of a dependence 
on the initial electron spectrum. We have experimented with a few 
different initial spectra to confirm this. 



We assume an isotropic distribution of electron pitch an- 
gles. To generate maps of projected synchrotron emis- 
sion, we construct a 2-D grid upon which the tracer par- 
ticle luminosities are mapped according to the "cloud- 
in-cell" (GIG) prescription (iHocknev fc Eastwood! [1988I ) 
and projected along the chosen line of sight. The result- 
ing synchrotron brightness for each sky pixel is given by 
Ij^ = (l + 2;)I/j^(i^2)/47rI)|^/Ar^, where Dl is the luminos- 
ity distance and Al] = /S.x/S.y/ D\^ where Da is the an- 
gular diameter distance. Our mock brightness maps are 
then convolved with a 2-D gaussian of FWHM 10 kpc, 
to simulate the effect of the PSF for a high-resolution 
instrument (this corresponds to 3" -10" for z = 0.05-0.2). 
No attempt was made to simulate the interferometric ef- 
fects on the images. 

Figures [10] through [12] show the projected gas tem- 
perature with radio brightness contours overlaid, for the 
frequencies 153, 327, and 1420 MHz in the 2:-projection 
(these fiducial values are selected to correspond to the 
GMRT and VLA frequencies). From these maps, it can 
be seen that the radio emission at lower frequencies per- 
sists over a long period and is bounded by the core cold 
fronts, apparent in the temperature maps. The mini- 
halo emission at u = 1.4 GHz becomes dimmer, cov- 
ers a smaller area on the sky, and becomes more patchy 
and amorphous. The time evolution will be discussed 
in more detail below. The maps in the z-projection are 
particularly striking when compared to the minihalo in 
RXJ 1720.1+26 from Figure [H 

In the projections that are in the orbital plane (along 
the X and ?/-axes), it is difficult to see the cold fronts 
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Fig. 11. — Projected gas temperature maps with radio contours overlaid at several epochs for the frequencies 153, 327, and 1420 MHz in the 
x-projection. The colorbar is temperature in keV. Contours of radio emission at (153, 327, 1420) MHz begin at (1.0, 0.5, 0.125) x 10""^ mJy 
arcsec""^ and increase by a factor of 2. Each panel is 750 kpc on a side. Tick marks indicate a distance of 100 kpc. 



in the temperature maps, but the radio emission is stih 
clearly bounded by them, with a radius r ~ 100 — 300 kpc 
that increases as the volume of the sloshing region in- 
creases with time. Figure [13] shows example profiles of 
the 327 MHz radio emission in the z-projection at the 
epoch t = 3.35 Gyr, demonstrating the lack of radio emis- 
sion beyond the cold front surfaces. The spatial coinci- 
dence of the radio emission with the X-ray cold fronts 
in the cool core, and the steep cutoff of the mini-halo 
emission are in agreement with observed mini-halos (e.g., 
[Mazzotta & Giacintucci 2008; Giacintucci et al. 2011). 

It is instructive to examine the properties of the ther- 
mal plasma over the same period to compare to the fea- 
tures seen in the radio maps. The radio emission will be 
dependent upon the turbulent velocity to continuously 
reaccelerate electrons, and the magnetic field to produce 
the emission itself. Figures [M] through [16] show the pro- 
jected magnetic field strength and turbulent velocity over 
the same epochs as the radio emission in the preced- 
ing figures. Within the volume of the sloshing region, 
the magnetic field has been significantly amplified, and 
the turbulence is strongest (though not necessarily in the 
same locations within the cool core). The fact that both 
of these effects are bounded within the envelopes of the 
cold fronts is what constrains the radio emission to these 
boundaries. Additionally, within the sloshing region, the 
amplification of the magnetic field and the strength of the 
turbulent velocity is far from uniform. The regions with 
the highest turbulent velocities {Sv ^ 100 — 200 km s~^) 
span spatial scales of ~50-100 kpc. Significant fluctua- 
tions in both of these quantities result in stronger emis- 



sion in localized regions. In particular, it appears that 
the regions of brightest radio emission correspond to the 
regions which have the highest turbulent velocities. This 
is due to the fact that the reacceleration timescale, which 
is on the order of 0.1 Gyr, is less than the timescale of 
the drift of regions with high turbulent motions, which 
is on the order of Gyr. 

One important observed characteristic of mini-halos is 
their steep radio spectra. This is naturally explained 
by the reacceleration model implemented in our simu- 
lations. The balance between reacceleration and losses 
sets a cutoff energy at which there is a sharp drop in the 
relativistic electron population, which in turn produces 
a steepening in the synchrotron spectrum. Since there 
is a range of magnetic field strengths in the cluster core, 
this break frequency will be different for different elec- 
trons, and the resulting spectrum will gradually steepen 
at higher frequencies. Assuming a power-law spectrum 
for the radio emission of the form Ij^ ex u~^^ all observed 
minihalos have a spectral index a ~ 1 — 2 around the 
frequency of 1 GHz, a much steeper spectrum than the 
radio emission from radio galaxies {a ~ 0.5 — O.S. ICondonl 
1992). 

Figure [17] shows the evolution of the synchrotron spec- 
trum of the emission within 300 kpc of the cluster center 
over the duration of the simulation. The initial spectrum 
at the epoch t = 2.55 Gyr is not shown on the plot since 
it is many orders of magnitude fainter than the spec- 
trum at later epochs. After the initial quick increase in 
emission on a 0.5-1 Gyr timescale, from the epoch t ~ 
3 Gyr onward the spectrum of the mini-halo becomes 
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Fig. 12. — Projected gas temperature maps with radio contours overlaid at several epochs for the frequencies 153, 327, and 1420 MHz in the 
y-projection. The colorbar is temperature in keV. Contours of radio emission at (153, 327, 1420) MHz begin at (1.0, 0.5, 0.125) x 10~^ mJy 
arcsec""^ and increase by a factor of 2. Each panel is 750 kpc on a side. Tick marks indicate a distance of 100 kpc. 
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Fig. 13. — Profiles of gas temperature and radio emission across cold front surfaces. Left: Projected gas temperature maps with radio 
contours overlaid at the epoch t = 3.35 Gyr, in the 2;-projection. The colorbar is temperature in keV. Contours are of 327 MHz radio 
emission which begin at 5 x 10~^ mJy arcsec"^ and increase by a factor of 2. The panel is 750 kpc on a side. Tick marks indicate 100 kpc 
distances. Right: Profiles of projected gas temperature and radio emission from the left panel, along the lines in the left panel. 



fainter and steepens at higher frequencies, as coohng be- 
gins to dominate over turbulence and fewer electrons are 
able to emit significant emission at these frequencies. To 
quantify this steepeing, we determine the spectral index 
of our simulated emission for these epcohs. We fit the 
total synchrotron spectrum to a power-law P^ ex u~^ 
over the frequencies 327, 610, and 1420 MHz, assum- 
ing the errors to be 10% of the fiux at each frequency. 



The evolution of the spectral index is tabulated in Fig- 
ure [171 For most of the evolution of the mini-halo, the 
spectral index hovers around a ^ 1 — 2, comparable to 
that of observed mini- hal o sources . The spectra for the 
minihalos in the Pe rseus (|Sijbring| [l993) and Ophiuchus 
(|Murgia et al.|[2QTQl ) clusters are shown for comparison. 
The spectrum of simulated cluster is compatible in both 
shape and normalization (though note that our normal- 
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Fig. 14. — Left: Projected (volume- weighted) magnetic field strength in the z-direction for the epochs t = 2.95, 3.35, 3.75, and 4.15 Gyr. 
Right: Projected (mass- weighted) turbulent velocity (estimated using only the Vz component and scaled to match the total turbulent 
energy) in the z-direction for the same epochs. Each panel is 750 kpc on a side. Major tick marks indicate 100 kpc distances. 





Fig. 15. — Left: Projected (volume-weighted) magnetic field strength in the x-direction for the epochs t = 2.95, 3.35, 3.75, and 4.15 Gyr. 
Right: Projected (mass- weighted) turbulent velocity (estimated using only the Vz component and scaled to match the total turbulent 
energy) in the x-direction for the same epochs. Each panel is 750 kpc on a side. Major tick marks indicate 100 kpc distances. 



ization is uncertain by an order of magnitude due to the 
weak constraints provided by IC measurements). 

Figure [TSl shows the evolution of the total synchrotron 
power within r = 300 kpc of the cluster center in the 
z-projection in W/Hz at the frequencies 153, 327, and 
1420 MHz. From the beginning of the injection of the 
relativistic particles at t = 2.55 Gyr until t ~ 3.2 Gyr, 
the synchrotron power increases from initially low values 
due to the acceleration of the low-energy particles. At 
t ^ 3.2 Gyr, the synchrotron power at each frequency 
reaches a maximum, in line with the evolution of the 
electron spectrum (see Figure [9]). The peak of the ra- 
dio power for our simulated cluster is consistent with 



that of observed mini-halos, which typically fall in the 
ran^e of Pi .4 qhz ^ a few xlO^^ - a few xlO^^ W/Hz 
(Cassa iiQ et al.|[2008L again, note that our normalization 



is uncertain by an order of magnitude.). After this, the 
power at each frequency slowly decays. The period dur- 
ing which the total radio power of the simulated mini- 
halo is within a factor of 10 from its peak value is about 
1 Gyr for 1.4 GHz, and slightly longer for lower fre- 
quencies (consistent with the evolution of the electron 
spectrum-the 1.4 GHz emission comes from electrons 
with 7 > (1 — 2) X 10^). This is significant, because it 
indicates that if mini-halos are powered by reacceleration 
they must be transient sources, particularly at higher fre- 
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Fig. 16. — Left: Projected (volume- weighted) magnetic field strength in the y-direction for the epochs t = 2.95, 3.35, 3.75, and 4.15 Gyr. 
Right: Projected (mass- weighted) turbulent velocity (estimated using only the Vz component and scaled to match the total turbulent 
energy) in the y-direction for the same epochs. Each panel is 750 kpc on a side. Major tick marks indicate 100 kpc distances. 
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Fig. 17. — Time dependence of the spectrum of the mini-halo 
within a radius of 300 kpc of the cluster center in W/Hz vs. the 
frequency of the emission in MHz. The effective power law slope 
a betwen the frequencies 327-1420 MHz is listed in the key, with 
la errors given. The arrows indicate the rise and fall of the spec- 
trum with tim e. The spectra f or the Perseus (Siibring 1993) and 
Ophiuchus (fMurgia et al.|[20Tol ) cluster minihalos are plotted for 
comparison. 



Fig. 18. — Time dependence of the synchrotron powers of the 
mini-halo in W/Hz at the frequencies 153, 327, and 1420 MHz 
within a radius of 300 kpc of the cluster center. The initial radio 
power at the frequency of 1420 MHz is much smaller than shown 
(P ~ 10-'^'-' W/Hz); the lower limit of the y-axis is set at a higher 
value for clarity. 



quencies. These properties are similar to those expected 
in the case of giant radio halos in the turbulent reac- 

celeration mod el. Based on GMRT radio survey data, 

iCassan o et alJ (^2008') concluded that mini-halos are rare, 
and interpreted this as support for their origin as a result 
of transient reacceleration events due to minor merging 
activity. Additionally, the longer lifetime of radio emis- 
sion at lower frequencies suggests that future surveys of 
clusters of galaxies at lower radio frequencies with obser- 
vatories such as LOFAR might reveal many more mini- 



halos in relaxed clusters. 



5.3. 



Simulated Inverse Compton Emission 

To illustrate the spatial distribution of relativistic elec- 
trons disentangled from the magnetic field distribution, 
we compute maps of the inverse- Compton (IC) emission 
for our tracer particles. The monochromatic IC power 
for a distribution of relativistic electrons as a function of 
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Fig. 19. — Maps of monochromatic IC surface brightness from the relativistic electrons in the simulation for two representative epochs 
for 50 keV photons. Each panel is 750 kpc on a side. Tick marks indicate 100 kpc distances. 



10"^ 
10^- 

10^- 



"'^'--^ll^^.. 


>^- 


— Thermal 
-- IC,t = 2.75 Gyr 
. -- IC, t = 2.95 Gyr 
-- IC,t = 3.35 Gyr 
IC,t = 3.75 Gyr 
-- IC,t = 4.15 Gyr 


---::§^ 







10 

E (keV) 




3.25 3.50 3.75 

t (Gyr) 



Fig. 20. — The monochromatic X-ray luminosity for several dif- 
ferent epochs of the simulation within a radius of 300 kpc from 
the cluster center. The black solid line indicates X-ray emission 
from the thermal gas at the epoch t = 2.55 Gyr, which has been 
computed using the APEC model ( Smith et al. 2001) assuming a 
spatially uniform metallicity of Z = 0.3 Zq. Dashed lines indicate 
IC emission, computed as described in the text. 



emitted photon energy ei is (|Rvbicki fc Lightmanlll979f ): 

P(ei) = Icarjcle (^) v{e) J djN{j)j-^f (^^ 

(37) 
where '^(e) is the incident photon number density at the 
incident photon energy e, ctt is the Thomson cross sec- 
tion, and 

f{x) = 2xlnx + x + l-2x^0 < x < 1. (38) 

Given that our distribution function for each tracer par- 
ticle is simply Nj{'y) = Kj X^i ^(7 ~ Tij)? ^^^ total IC 
power at a given frequency for each tracer particle is 



Fig. 21. — The evolution of hard X-ray luminosity integrated 
over the 20-80 keV waveband for the emission within a radius of 
300 kpc from the cluster center. 



then given by: 



Pj{ei) = 3caTK^ 
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(^) >: h^ / h^ de. 



H.^ 






(39) 

We assume that v{e) is a blackbody spectrum with tem- 
perature ksTcMB (which is redshift-dependent) and in- 
tegrate over all photon energies and sum over all jij to 
obtain the emitted power at each photon energy. 

Since our main aim in this paper is to determine the 
properties of the radio emission of the minihalo, our ex- 
amination of the properties of the IC emission will be 
comparatively brief. Figure [19] shows maps of the IC 
intensity at the epochs t = 2.95 and 3.75 Gyr at the 
representative photon energy E^ = 50 keV. As in the 
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case of the radio emission, the hard X-ray emission at 
high energies is bounded by the cold front surfaces. The 
spatial distribution of the IC emission is similar to the 
synchrotron emission. This is expected, since, as we men- 
tioned in the previous section, the spatial distribution of 
the radio emission follows the location of the turbulent 
regions, and both kinds of emission originate from the 
same populations of relativistic electrons. 

Figure [20] shows the IC spectrum at a few different 
epochs of the simulation for the central 300 kpc of the 
cluster core, compared to the thermal spectrum at the 
epoch t = 2.55 Gyr from the same core region (a mix- 
ture of temperatures in the range T ~2-5 keV). At nearly 
all energies E^ ~ 1 — 100 keV, the spectrum is dominated 
by the thermal emission. It will not be possible to de- 
tect such emission against the thermal emission from the 
bright cluster core with upcoming X-ray telescopes with 
hard X-ray detection capabilities such as NuSTAR and 
ASTRO-H^ except possibly at energies near 100 keV. 

For synchrotron and IC emission produced by the 
same population of relativistic electrons, Psync/Pic = 
^^/^CMB- Taking the volume- averaged magnetic field 
By ^ 5 /iC and Bqmb ^ 4 fiG (assuming z = 0.1), 
we find Psync/Pic ^ 1-56. Confirming this, we estimate 
from Figures [TTl and [2Ql that the two bolometric luminosi- 
ties are roughly equal with Pgync ^ Pic ^ 10^^ erg s~^, 
which is orders of magnitude smaller than the ther- 
mal X-ray luminosity of Ptherm ^ 10^^ erg s~^. This 
would be the case for all clusters hosting a minihalo, 
as most cool-core clusters have relati vely s trong mag- 
netic fields in their cores (iTavlor etS]|2QQ2i I2006L 120071 : 
iBonaf ede et al.l l2010l : IVacca et al.ll2011f ). 

Figure [2T] shows the evolution of the X-ray luminosity 
in the energy band 20-80 keV, for the r < 300 kpc region 
of the cluster core. The behavior of the hard X-ray power 
is similar to that of the synchrotron power in that it ex- 
hibits a fast increase (of nearly an order of magnitude) 
when reacceleration begins, followed by a slower dropoff. 
The peak of the hard X-ray power is ~ 2 x 10^^ erg s~^, a 
few orders of magnitude smaller than the non-thermal lu- 
minosity possibly detected by XMM-Newt on and INTE- 
GRAL i n Ophiuchus ^ 10^ ^ erg s"^ (^Murgia et al.ll2010L 
also see lAiello et all (120091 ) for upper limits from Swift 
and iFujita et al.l (|2008[ ) for upper limits from Suzaku) . 



5.4. Different Initial Spatial Distributions of 
Relativistic Electrons 

Though our tracer particles are distributed through- 
out the simulation domain, it is not necessarily the case 
that in real clusters relativistic electrons are similarly dis- 
tributed. Electrons are accelerated to relativistic speeds 
by AGN and generated as byproducts of collisions of rela- 
tivistic protons with the thermal protons of the ICM. The 
former process will result in relativistic electrons prefer- 
entially located in AGN-blown bubbles, and the latter 
will produce more relativistic electrons in regions of high 
gas density throughout the cluster. The appearance of 
a mini-halo generated by sloshing may depend on the 
initial spatial distribution of these particles. 

Our already examined case of a large spherical dis- 
tribution of relativistic particles represents the second 
possibility mentioned above. In order to approximately 
represent a population of seed electrons produced by an 
AGN, we perform three simulations of relativistic elec- 



trons originating in two "bubbles". In these cases, it 
order to produce radio emission that follows the X-ray 
cold fronts, it will not only be necessary for electrons to 
be reaccelerated by the turbulence, but the sloshing mo- 
tions will have to adequately redistribute the seed elec- 
trons throughout this region. 

Since modeling the dynamics of the relativistic AGN 
bubbles is far beyond our current scope, we perform a 
simple exercise that tests whether or not the relativistic 
particles may be adequately redistributed by the sloshing 
motions to fill the core region. We first assume that the 
relativistic particles contained in the bubbles have just 
become mixed in with the thermal gas at the epoch t = 
1.95 Gyr (before the beginning of the sloshing motions), 
and then identify tracer particles contained within two 
symmetric 30 kpc "bubbles" centered around the clus- 
ter potential minimum (without attempting to evacuate 
thermal gas from these "bubbles" , as would be the case 
in real clusters). Each of these particles are assumed to 
contain the same mass of gas rrij as before, but there 
are far fewer of them than in our default spatial case, 
with approximately ~ 10^ tracer particles per bubble. 
We then follow these particles for 0.6 Gyr, and at t = 
2.55 Gyr, we assign them relativistic particle distribu- 
tions identical to that of our default case. We perform 
three simulations of the evolution of the relativistic parti- 
cles from these initial distributions, with the axis of the 
bubbles aligned along the x, y^ and z axes (keeping in 
mind that in our simulation, sloshing occurs in the x — y 
plane. 

Figure [22] shows the evolution of particle positions be- 
tween the epochs t = 1.95 Gyr and 2.55 Gyr for our 
three bubble initial conditions. Regardless of the initial 
orientation of the bubbles, we find that the sloshing mo- 
tions redistribute the gas particles so that most of them 
end up within the spiral shape traced out by the cold 
fronts. Figure [23] shows the projected gas temperature 
in the z-direction of the simulation with 327 MHz ra- 
dio contours overlaid at the epoch t = 2.95 Gyr for two 
simulations where "bubbles" of relativistic particles have 
been injected. We find that faint radio emission is pro- 
duced in all three cases, and it is contained within the 
core of the sloshing region. Although the particles have 
been redistributed throughout the sloshing region, only a 
fraction of them have encountered regions of strong tur- 
bulence, fewer than in the case where we assumed rela- 
tivistic particles were originally distributed proportional 
to gas density everywhere within the core. 

If we assume that a population of low-7 electrons is 
built up over time in the core by injection from the cen- 
tral AGN and then mixing by sloshing motions, then we 
may end up with a situation similar to our default case 
for the spatial distribution of relativistic electrons, which 
succeeded in reproducing the observed mini-halo proper- 
ties. This lends support to the hypothesis that mini- 
halos may be generated by electrons which originated 
from the central AGN and are reaccelerated by sloshing 
motions. Sloshing and AGN are both commonly found 
in cool-core clusters where mini-halos are found. Addi- 
tionally, as po inted out by Pfrommer & EnBlin (20( 3); 
ICassanoet aH (|2008[) : FKeshet fc Loeh (.2010. ): .Keshetl 
(I2OIOI ), in the cool cores of galaxy clusters relativistic 
electrons will be efficiently produced by hadronic pro- 
cesses due to the high density of thermal protons. Be- 
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Fig. 22. — Evolution of the positions of the gas particles between the epochs t = 1.95 Gyr and 2.55 Gyr for the three different initial 
spatial configuations of "bubble" particles, overlaid on maps of projected gas temperature. Each panel is 200 kpc on a side. Major tick 
marks indicate 50 kpc distances. 




Fig. 23. — Projected gas temperature maps with radio contours overlaid from different "bubble" initial conditions for the relativist ic 
particles. The colorbar is temperature in keV. Contours are of 327 MHz radio emission which begin at 1.25 x lO"'^ mJy arcsec""^ and 
increase by a factor of 2. Left: initial bubbles along x-axis. Center: initial bubbles along y-axis. Right: initial bubbles along z-axis. The 
epoch is t = 2.95 Gyr. Each panel is 400 kpc on a side. Tick marks indicate 100 kpc distances. 



tween the AGN and hadronic interactions, the dense 
cores of clusters should contain adequate numbers of seed 
relativistic electrons to be reaccelerated by turbulence. 

6. DISCUSSION 

6.1. The Evolution of the Relativistic and Turbulent 
Energy Components 

It is instructive to examine the energy budget of the 
various components of our model as they evolve during 



the simulation. The solid lines in Figure [24] shows the 
evolution of the total kinetic, turbulent kinetic, mag- 
netic, and relativistic electron energies relative to the to- 
tal thermal energy of the tracer particles initially within 
a radius of 100 kpc from the cluster center. The largest 
contribution is that of the total kinetic energy, which 
stays around ^3-4% over time, followed by the mag- 
netic energy at ~1% and the turbulent kinetic energy 
at ~0.5%. The relativistic electron energy density starts 
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Fig. 24. — Time dependence of the various energy components 
of the simulation. Sohd hues indicate quantities computed for all 
of the tracer particles originating within a radius of r = 100 kpc at 
the epoch f = 2.55 Gyr. Dashed lines indicate quantities computed 
for all of the tracer particles which at some point in the simulation 
fulfill the criterion 7 > 10^. 



at our initial value of £rei/^th — 10~^ and decreases as 
most of the electrons cool over time-the liigli-7 tail of the 
electron distribution that is studied here does not con- 
tribute much to the total electron energy. All of the sep- 
arate nonthermal energy categories combined contribute 
less than 10% to the total energy budget. This is in 
part a consequence of the specialized initial conditions 
of our setup, which began with only thermal, magnetic, 
and a small amount of kinetic energy (due to our relax- 
ation procedure, see Section [2TT]) . This is consistent wit h 
the cosmological simulation results of iLau et al.l (|2QQ9l ). 
who show that for cores of relaxed clusters, the contribu- 
tion to the total pressure from turbulence, rotation, and 
streaming gas motions is negligible compared to the ther- 
mal pressure, as well as recent constraints placed on tur- 
bulent motions in the core of Abell 1835 bylSanders et aP 
(|2QTQh . and in Abeh 3112 bv lBulbul et"atl ([Ml ). 

It is also interesting to examine the energies of the 
tracer particles that are emitting the most synchrotron 
radiation. These particles will be the ones that have en- 
countered the highest turbulent velocities, as can be seen 
by the correspondences between the regions of brightest 
synchrotron emission (Figures [TQ]fT2]) and the locations 
with the strongest turbulence (Figures [T4lfT5]) . To isolate 
these regions, we select only those tracer particles that 
will reach an average energy 7 > 10^ at any point in 
the simulation and emit at frequencies ^ 100 MHz, at 
which the halos are presently observed. These particles 
comprise 1% of the total number of tracer particles of our 
initial set; their energy ratios are shown by the dashed 
lines of Figure [24l The synchrotron emission indeed high- 
lights regions with slightly higher turbulent energies and 
(by virtue of their selection) nearly an order of magni- 
tude higher relativistic electron energy contribution. 

We can also determine the energy gained by relativis- 
tic particles through turbulent acceleration. Though we 
have not modeled self-consistently the interaction be- 



tween the relativistic particles and the turbulent ICM 
(in particular their back-reaction on turbulence), for the 
reacceleration model to be a viable mechanism for gener- 
ating mini-halos, the energy gain by relativistic particles 
should be a small fraction of the turbulent energy. To 
do this, we need to compare the rate of energy gained by 
the relativistic electrons versus the rate of energy that 
passes through the turbulent cascade toward the dissi- 
pation scale. The former rate will not be equal to the 
total change in energy in relativistic particles, as most of 
the energy gained (especially for the highest-energy elec- 
trons) will end up in the form of synchrotron and hard X- 
ray emission, but we know from Equation [18] how much 
energy the electrons in our simulations gain. To get a 
rough determination of the ratio of these two quantities, 
we compute the former via Equation [T71 and the latter 
via the estimate f or the turbul ent energy cascade rate 
(Equation 48 from iBrunetti fcXazarian 2007) : 



^turb ^ C 



(M)^(c. 



(40) 



where ^ = 30 kpc and the constant C ~ 5 — 6. In making 
this estimate we have used the average turbulent Mach 
number M. — Vt/cg and sound speed Cg. The ratio of 
these rates for all the tracer particles and the high-7 
tracer particles is also shown in Figure [231 The rate of 
energy gain of the relativistic particles is typically ~10~^ 
of the turbulent cascade rate. 

As a consistency check, we check the ratio of the rate 
of the turbulent energy cascade rate to the cooling rate 
of the thermal gas by X-rays. We find that the rate of 
turbulent dissipation is at the level of ~10% of the bolo- 
metric X-ray luminosity, however, this estimate is prob- 
ably uncertain by an order of magnitude, since we have 
made an approximate estimate of the turbulent cascade 
rate due to our inability to measure this quantity from 
the simulation directly. Thus, this quantity is probably 
overestimated; a consequence of which is that the effi- 
ciency of reacceleration relative to the turbulent cascade 
rate is likely to be somewhat higher than this estimate 
indicates. 

In general, only a fraction of the energy of turbu- 
lence is drained int o relativistic particles via t urbulent 
reacceleration (see IBrunetti fc Lazarian] l2011b[ ). Con- 
sequently, even if we have not considered the case of 
cosmic ray protons that could be present in the clus- 
ter core region, the additional energy input that would 
be acquired by these pre-existing protons via accelera- 
tion due to sloshing-driven turbulence in our simulation 
is ^ 1% of the thermal energy, regardless of the initial 
energy content of the pre-existing protons. A energy 
content of cosmic ray protons at percent level of the 
thermal gas is still consistent with pre sent upper limits 
from gamma-ray observations (see, e.g. Aharonian e t al.l 
2009; Ackermann et al. 2010; Jelt e ma fc Profum o 20TTI : 
IMAGIC Collaboration efliIll20TTI) . 

6.2. Two Effects of Sloshing on the Relativistic 
Electrons 

From our simulations, we discern two effects of the gas 
sloshing on the relativistic electrons that are important 
for the formation of mini-halos. The first is the reac- 
celeration itself. Without reacceleration, the electrons 
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will simply cool down to energies below which they are 
incapable of emitting at observable synchrotron frequen- 
cies. With reacceleration, the electrons inside the slosh- 
ing region are maintained for a significant period of time 
(~1 Gyr from the onset of sloshing) at high enough 7 
for radio emission to be produced (see Figure [9]). The 
second important effect of sloshing concerns the initial 
spatial distribution of the electrons. Sloshing motions 
can take rather general initial spatial distributions of 
relativistic electrons and redistribute them throughout 
the sloshing region. Thus, we anticipate that for reason- 
able seed populations of relativistic electrons (either re- 
sulting from hadronic interactions or injected by AGN), 
sloshing should be able to redistribute and reaccelerate 
them so that radio mini-halo emission is produced. As 
previously mentioned, sloshing also amplifies magnetic 
fields within the cluster core, an effect that is als o im- 
portant for hadron ic /secondary models (Keshet fc Loebl 
l2(M lKeshetl l20T0h. 

6.3. Morphology of Mini-Halos 

Some mini-halos in observed galaxy clusters with 
attendant sloshing cold fronts have a very specific 
morphology-they are very similar to the shape of the 
sloshing region and bounded by the cold fronts as seen in 
X-rays (Mazz otta fc Giacintucci 2008) . However, other 
mini-halos do not appear to be coincident with sloshin g 
cold fronts (.Govoni et al.l 120091 : [Giacintucci et al.ll2011[ ). 
and some of them appear in clusters without discernable 
sloshing. However, turbulence generated by sloshing may 
still be the origin of the radio emitting electrons, even if 
the corresponding cold fronts are not seen. 

Projection effects may make it difficult to associate 
mini-halos with cold fronts. In projections perpendic- 
ular to the plane of the mutual orbit of the main cluster 
and the subcluster, such as the z-projection of our sim- 
ulations, a clear association of the radio emission with 
the shape of the cold front is observed. In other pro- 
jections {x and ?/), the spiral shape of the cold fronts is 
unobservable, though small portions of the cold fronts 
still appear. In these cases, the emission is still bounded 
within this region, but it is more difficult to tell from the 
X-ray image if sloshing is present. 

Larger cold fronts will also be more difficult to observe 
in X-rays due to their low surface brightness contrast. 
Finally, although the conditions necessary to generate 
the mini-halo emission may be associated with the slosh- 
ing cold fronts, it does not necessarily imply that these 
conditions prevail uniformly throughout the sloshing re- 
gion. Figures [14] through [T6] show the turbulent kinetic 
energy exhibits strong variations throughout the cluster 
core. Though the magnetic field strength is the strongest 
along the cold front surfaces, it also varies greatly inside 
them, also seen in Figures [14] through [16] These varia- 
tions may combine to produce emission maps that do not 
always fill the entire region bounded by the cold fronts 
(see Figures [To] through [T2]). particularly since the spatial 
scale of these variations is comparable to the size of the 
core region. 

6.4. Uncertainties in the Model 

Finally, it is important to summarize the various un- 
certainties involved in our model. The major sources 



of uncertainty arise from the assumptions of the initial 
energy in relativistic electrons and the assumptions in- 
volved in the estimates of the reacceleration coefficients. 

Our initial electron spectrum depends on a number of 
conditions as detailed in Section 13.11 We have experi- 
mented with a few different initial electron spectra un- 
der otherwise identical simulation conditions, and have 
found that our results do not depend significantly on the 
initial spectral shape. All of our results related to the rel- 
ativistic electron spectra, including our synchrotron and 
IC emission, may be scaled by the relativistic electron en- 
ergy fraction rj. The hard X-ray constraints o n the clus- 
ter IC emission allow values up to r] ~ 10~^ (|Wik et al.l 
I2OI2I ): a value in equipartition with the magnetic field 
energy would also be 77 ~ 10~^; thus our use of 77 = 10~^ 
is very conservative. 

The reacceleration coefficients I^pp,TTD and ^pp,c in 
our model are dependent upon the characteristics of the 
power spectrum of turbulence, in particular the mini- 
mum and maximum wavenumbers /cmin and /cmax, and 
the ratio of power in compressible motions to the total 
power R^. While kmin is determined by the size of the 
sloshing structures and cannot be much different from 
our assumption, the value of kcut represents an exten- 
sion (albeit a very plausible and likely conservative one) 
by a factor of ~10 above the range of wavenumbers for 
the inertial regime in our simulations. If our kcut is an 
overestimate, we subsequently overestimate Dpp and the 
reacceleration coefficient x by a factor of '^ a few. Con- 
versely, if the turbulent cascade extends down to sni aller 
scales in the manner of lBrunetti fc LazarianI (|2011bl ) , we 
have underestimated kcut by a factor of ^ a few, and the 
reacceleration efficiency would be substantially underes- 
timated. 

We already discussed the potential overestimate in Re 
due to the edge effects of the non-periodic subdomains 
in Section 14. 2[ which showed that up to ~30% of the 
compressive power is potentially spurious. However, the 
main uncertainty in i?^ arises from our adoption of an 
average value constant in space and time. In reality, the 
compressive component of turbulence will vary spatially 
and with time. Our adoption of R^ = 0.25 is a conserva- 
tive one, since most of the reacceleration from turbulence 
will arise from regions with higher i?^, which we found 
can be as high as 0.7 (Section l4.2p . In these regions, the 
reacceleration due to turbulence is underestimated by a 
factor of a few. 

Figure [25] shows the effects of varying the coefficient 
Dpp by a small factor in either direction on the total 
electron spectrum of the simulation. Decreasing Dpp by 
a factor of two results in a decrease of A/" (7) by an order of 
magnitude or more at 7 ^ 10^. Increasing Dpp by a factor 
of 1.5 results in an increase of N{j) by nearly an order of 
magnitude or more at 7^10^. This demonstrates that 
the resulting electron spectrum is fairly sensitive to the 
value of Dpp at high 7. 

Due to the combined effects of these uncertainties, our 
procedure must be regarded only as an order of magni- 
tude estimate, which is suitable for the aim of the present 
paper. Nevertheless, for the physically plausible choices 
we have made for these coefficients, our results are con- 
sistent with observations of minihalos. 
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Fig. 25. — Dependence of the total electron spectrum on the value 
of Dpp, with values reported relative to the default value adopted 
in the paper. 



7. SUMMARY 

In this work, we have performed the first simulation 
test of the hypothesis that radio minihalos, sometimes 
observed in cluster cool cores, originate from relativistic 
electrons reaccelerated by turbulence generated by slosh- 
ing motions. The seed low-7 electrons are likely to re- 
main, e.g., from past AGN activity and/or from hadronic 
interactions. The acceleration mechanisms that we con- 
sidered are the damping of turbulence-induced magne- 
tosonic waves on relativistic electrons and non-resonant 
compression, which are included in a subgrid fashion. 
Our high-resolution, idealized simulation consists of a 
cool-core galaxy cluster with magnetized gas and an in- 
falling gasless subcluster. The interaction between the 
subcluster and the core of the main cluster initiates gas 
sloshing in the core. Two significant magnetohydrody- 
namic effects of the gas sloshing are the generation of 
turbulence and the amplification of magnetic fields, both 
of which occur within the region of the cool core bounded 
by the sloshing cold fronts. 

We have shown that in this model, faint, extended ra- 
dio sources that closely resemble the spatial and spec- 
tral properties of the observed mini-halos may be gener- 
ated under a variety of physically reasonable assumptions 
for the initial spatial distribution and input spectrum of 
seed relativistic electrons. The sloshing motions-bulk 
fiows and turbulence-have two effects on the relativistic 
electrons, reacceleration and redistribution of seed rel- 
ativistic electrons throughout the core. These two ef- 
fects combine to produce radio emission that is diffuse, 
steep-spectrum, and that traces the spatial features of 
the X-ray emitting gas. The radio power generated in 
our model is in the range of that of observed miniha- 
los (though our assumption of the initial energy density 



of the relativistic electrons is rather uncertain). Impor- 
tantly, we have shown that minihalos produced by tur- 
bulent reacceleration are transient sources, particularly 
at high (~1 GHz) frequencies, consistent with the small 
number of observed mini-halos. We have also shown 
that the X-ray inverse- Compton emission produced by 
the same relativistic electrons will be difficult to observe 
due to the fact that the IC emission begins to dominate 
over the thermal emission of the bright core only at very 
high energies (~100 keV). 

Significant improvements could be made to bring our 
model more in line with observations. More sophisti- 
cated methods can be developed to better characterize 
the spectrum and normalization of turbul ence created 
by sl oshing (such as the recent attempts in IVazza et al.l 
I2QI2I ). A more accurate treatment of the evolution of 
the spectrum of relativistic electrons requires the use of 
a physical model for a time-dependent injection of these 
electrons. Also, in our paper we do not consider the 
effect of cosmic ray protons in the sloshing region. Self- 
consistent calculations of particle acceleration and of the 
resulting non-thermal emission should take into account 
also the proton component, since protons generate sec- 
ondary electrons (that eventually may be reaccelerated) 
and can be important for the damping of the turbulence 
(jBrunetti fc Lazarianll2Qlla[ ). Finally, a variety of clus- 
ter initial conditions could be explored, to determine if 
mini-halo formation is more likely under certain merger 
conditions than others. We leave these considerations for 
future work. Our simulation data could also be used to 
examine the hadronic hypothesis for mini-halos, and pro- 
vide a point of comparison between the two models. We 
will make this comparison the subject of a future paper. 

Our present work suggests that reacceleration of rela- 
tivistic electrons by turbulence can produce radio mini- 
halo emission in cool-core clusters, though the uncertain- 
ties at various steps of our simulations are high. 
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Fig. 26. — Verification test of simulation of relativist ic electrons with losses only. Parameters are B 
: = 0.1. Blue curve: input electron spectrum. Red curve: predicted electron spectrum. 
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APPENDIX 
TEST CASES FOR THE RELATIVISTIC INTEGRATOR 

In this appendix we present the results of a few verification tests of our relativist ic particle integrator. We assume in 
these test cases time-independent conditions for the state of the gas and the magnetic field. Additionally, for timescales 
that are short compared to the Hubble time, the effect of the redshift dependence on the CMB energy density is small 
and it may be ignored. We will first perform tests of the deterministic "drift" term of Equation [231 against analytic 
solutions, and follow with comparisons of a full solution of Equation [23] with the stochastic term included with a 
solution generated under identical conditions by integrating a Fokker- Planc k equation. 

Our tests of the deterministic evolution of particles are similar to those of .Sarazid (|1999l ). who derived the spectrum 
of cosmic-ray electrons in clusters of galaxies under the effects of synchrotron, IC, and Coulomb losses. For time- 
independent conditions the particle energy spectrum N{j^t) may be derived from an initial spectrum N{ji^ti) by a 
simple relation: 

N{j,t) = N[j,{rt),U]^-^^j^ (Al) 

where 6(7) is the rate of change of particle energy and 7i(7, i) is given impHcitly by 



W) 



{t - U) 



(A2) 



For our tests, we adopt parameters similar to those in [Sarazid ()1999l ). We evolve an initial power-law spectrum 
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Fig. 27. — Verification test of simulation of relativistic electrons with reacceleration. Parameters are 5v = 200 km s"-*^, B = 1 /xG, nth 
= 10~^, evolved from z = 0.1. Blue curve: input electron spectrum. Red curve: predicted electron spectrum. Green curve: predicted 
spectrum without reacceleration. The predicted "break frequency" 75 is shown by the dashed line. 
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Fig. 28. — Simulation of relativistic electrons with losses only for differing numbers of samples A^rel- Parameters are B = 1 /iG, nth = 
10~^, evolved from z = 0.1. The red curve indicates the predicted electron spectrum. Left panel: A^rel = 10^- Center panel: A^rel = 10^. 
Right panel: A^rel = 10^. 
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Fig. 29. — Simulation of relativistic electrons using the SDE integrator compared with solutions from a Fokker-Planck simulation. 
Parameters are B = Bcmb = 3.25 jiG, nth = 10""^, and tacc = X~'^ = 0-3 Gyr. 



(A^(7) oc 7~^) with p = 2.3 from z = 0.1 with B = 1 /iG, ng = 10~^ cm~^, and a cosmology with h = 0.65, ^rn = 1, 
and (Ia = 0. The evolved particle distributions are binned into energy bins of A7 = 10 to derive a particle energy 
spectrum N{'y) that may be compared with analytical results derived from Equations I AHIX2I Unless otherwise noted, 
the simulations used 10^ samples for the distribution function. 

We first consider the evolution of an initial population of relativistic electrons with no reacceleration. Figure [26l 
shows the resulting simulated electron spectrum A^(7) at 2; = 0, compared to the true spectrum (red line) and the 
initial spectrum at z = 0.1 (blue line). This spectrum is very similar to the electron spectra in Figures 6-8 of Sarazinj 
([1999 ). The effect of the Coulomb losses is to flatten the relativistic electron spectrum at low 7, whereas the effect of 
the synchrotron and IC losses is to introduce a sharp cutoff at high 7. 

Secondly, we include the effects of reacceleration. For this we use Equation [15] and assume a turbulent speed 
Svt = 200 km s~^ at the scale of irnin = 100 kpc. This spectrum is shown in Figure [271 compared to the true 
spectrum (red line), the initial spectrum (blue line), and the spectrum without reacceleration (green line). The 
resulting spectrum has a similar shape to the case without reacceleration, but it has been shifted to higher energies 
and the normalization of the spectrum is higher. If we assume reacceleration and losses are balanced at high 7, we 
may derive a "break frequency" 75 = x/f^ from Equations [T9l[2Ql For these conditions, 75 ^ 3131, which is reproduced 
well by the simulated data, as seen in Figure [271 

We also test the effects of varying the number of particle samples of our distribution function on the evolved spectrum. 
Figure [28l shows the evolved particle spectrum for three different values of A/'rei, assuming only particle energy losses 
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Fig. 30. — Simulation of relativistic electrons using the SDE integrator with the stochastic term included compared with a simulation 
where only the systematic terms are modeled. Parameters are B = Bqmb = 3.25 /iG, nth = 10""^, and tacc = X = 0-3 Gyr. 



and no reacceleration. For N^qi = 10^, the shape of the particle spectrum is barely discernable. For A/'rei = 10^, the 
shape is more apparent and for A^rei = 10^ the model is well-described by the simulated data. For the simulations 
presented in this work, we have chosen A^rei = 10^ as the number of samples per tracer particle as a good balance 
between accuracy and computation time. 

Finally, we verify our integration of the full stochastic differential equation by comparing the resulting electron 
energy spectrum to that produced by a Fokker-Planck calculation under identical conditions. For this test, we have 
assumed a constant magnetic field B = Bqub = 3.25 /iG, density of thermal particles nth = 10~^, and a reacceleration 



timescale of ta 



X" 



0.3 Gyr. As in the previous tests, the value of ^cmb was held fixed. The initial condition 



for the relativ istic particle distribution is the steady-state condition with continuous injection and energy losses from 
iSarazinI (|1999'), under the above conditions and a spectral index for the input spectrum of p = 2.5. 

Figure [29] shows the electron spectra produced by the two methods at the epochs of t = 0.16 and 0.48 Gyr. The 
agreement between the two methods over a wide range of 7 is excellent, except where 7 ^ 100 (where our method 
underpredicts A^(7) by ~10-40% depending on the epoch. In any case, this does not affect our results since the energies 
of these particles are too small to radiate at the radio frequencies of interest in this paper. 

To demonstrate the importance of including the stochastic term in our SDE, Figure [30] shows the evolved electron 
spectrum at the epoch t = 0.48 Gyr in a simulation where the stochastic effects are included with one where they are 
not, with identical conditions as in Figure [29] The effect of including the stochastic term in our model is significant 
for our results, as the number of electrons at the energies (7 > 10^) required for minihalo emission at the observed 
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frequencies is higher by orders of magnitude. 
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